Longitudinal monitoring of neutral and adaptive genomic diversity in a reintroduction.

Restoration programs in the form of ex-situ breeding combined with reintroductions are becoming critical to counteract demographic declines and species losses. Such programs are increasingly using genetic management to improve conservation outcomes. However, the lack of long-term monitoring of genetic indicators following reintroduction prevents assessments of the trajectory and persistence of reintroduced populations. We carried out an extensive monitoring program in the wild for a threatened small-bodied fish (southern pygmy perch, Nannoperca australis) to assess the long-term genomic effects of its captive breeding and reintroduction. The species was rescued prior to its extirpation from the terminal lakes of Australia's Murray-Darling Basin, and then used for genetically informed captive breeding and reintroductions. Subsequent annual or biannual monitoring of abundance, fitness, and occupancy over a period of 11 years, combined with postreintroduction genetic sampling, revealed survival and recruitment of reintroduced fish. Genomic analyses based on data from the original wild rescued, captive born, and reintroduced cohorts revealed low inbreeding and strong maintenance of neutral and candidate adaptive genomic diversity across multiple generations. An increasing trend in the effective population size of the reintroduced population was consistent with field monitoring data in demonstrating successful re-establishment of the species. This provides a rare empirical example that the adaptive potential of a locally extinct population can be maintained during genetically informed ex-situ conservation breeding and reintroduction into the wild. Strategies to improve biodiversity restoration via ex-situ conservation should include genetic-based captive breeding and longitudinal monitoring of standing genomic variation in reintroduced populations.

Monitoreo Longitudinal de la Diversidad Genómica Neutral y Adaptativa en una Reintroducción Marshall et al. 21-643 Resumen Los programas de restauración a manera de reproducción ex situ combinada con reintroducciones se están volviendo críticos para contrarrestar las declinaciones demográficas y la pérdida de especies. Dichos programas usan cada vez más la gestión genética para mejorar los resultados de conservación. Sin embargo, la falta de monitoreo a largo plazo de los indicadores genéticos posteriores a la reintroducción evita que se realicen evaluaciones de la trayectoria y la persistencia de las poblaciones reintroducidas. Se rescató un pez de talla pequeña (percha pigmea del sur [Nannoperca australis]) previo a su extirpación de los lagos terminales de la Cuenca Murray-Darling en Australia para después reproducirlo en cautiverio con información genética y reintroducirlo. Realizamos monitoreos anuales o bianuales de la abundancia, aptitud y ocupación en vida silvestre durante once años, además de un muestreo genético posterior a la reintroducción. Analizamos los datos genómicos de los grupos originales rescatados, los nacidos en cautiverio y los reintroducidos. Nuestro objetivo era evaluar los efectos genómicos a largo plazo de la reproducción en cautiverio y la reintroducción de esta especie. Esto reveló baja endogamia y el sólido mantenimiento de la diversidad genómica neutral y adaptativa durante varias generaciones. Encontramos una coherencia entre la tendencia creciente en el tamaño de la población efectiva de la población reintroducida y los datos de campo que demostraron el restablecimiento exitoso de la especie. Nuestro estudio proporciona un raro ejemplo empírico de cómo el potencial adaptativo de una población localmente extinta puede mantenerse durante la reproducción de conservación ex situ genéticamente informada y su reintroducción. Las estrategias para mejorar la restauración de la biodiversidad por medio de la conservación ex situ deberían incluir la reproducción en cautiverio basada en la genética y el monitoreo longitudinal de la variación genómica actual de las poblaciones reintroducidas.

Fish out of water: Genomic insights into persistence of rainbowfish populations in the desert.

How populations of aquatic fauna persist in extreme desert environments is an enigma. Individuals often breed and disperse during favorable conditions. Theory predicts that adaptive capacity should be low in small populations, such as in desert fishes. We integrated satellite-derived surface water data and population genomic diversity from 20,294 single-nucleotide polymorphisms across 344 individuals to understand metapopulation persistence of the desert rainbowfish (Melanotaenia splendida tatei) in central Australia. Desert rainbowfish showed very small effective population sizes, especially at peripheral populations, and low connectivity between river catchments. Yet, there was no evidence of population-level inbreeding and a signal of possible adaptive divergence associated with aridity was detected. Candidate genes for local adaptation included functions related to environmental cues and stressful conditions. Eco-evolutionary modeling showed that positive selection in refugial subpopulations combined with connectivity during flood periods can enable retention of adaptive diversity. Our study suggests that adaptive variation can be maintained in small populations and integrate with neutral metapopulation processes to allow persistence in the desert.

Movements and behaviour of blue whales satellite tagged in an Australian upwelling system.

Knowledge about the movement ecology of endangered species is needed to identify biologically important areas and the spatio-temporal scale of potential human impacts on species. Blue whales (Balaenoptera musculus) are endangered due to twentieth century whaling and currently threatened by human activities. In Australia, they feed in the Great Southern Australian Coastal Upwelling System (GSACUS) during the austral summer. We investigate their movements, occupancy, behaviour, and environmental drivers to inform conservation management. Thirteen whales were satellite tagged, biopsy sampled and photo-identified in 2015. All were genetically confirmed to be of the pygmy subspecies (B. m. brevicauda). In the GSACUS, whales spent most of their time over the continental shelf and likely foraging in association with several seascape variables (sea surface temperature variability, depth, wind speed, sea surface height anomaly, and chlorophyll a). When whales left the region, they migrated west and then north along the Australian coast until they reached West Timor and Indonesia, where their movements indicated breeding or foraging behaviour. These results highlight the importance of the GSACUS as a foraging ground for pygmy blue whales inhabiting the eastern Indian Ocean and indicate the whales' migratory route to proposed breeding grounds off Indonesia. Information about the spatio-temporal scale of potential human impacts can now be used to protect this little-known subspecies of blue whale.

Genome-wide association study of an unusual dolphin mortality event reveals candidate genes for susceptibility and resistance to cetacean morbillivirus.

Infectious diseases are significant demographic and evolutionary drivers of populations, but studies about the genetic basis of disease resistance and susceptibility are scarce in wildlife populations. Cetacean morbillivirus (CeMV) is a highly contagious disease that is increasing in both geographic distribution and incidence, causing unusual mortality events (UME) and killing tens of thousands of individuals across multiple cetacean species worldwide since the late 1980s. The largest CeMV outbreak in the Southern Hemisphere reported to date occurred in Australia in 2013, where it was a major factor in a UME, killing mainly young Indo-Pacific bottlenose dolphins (Tursiops aduncus). Using cases (nonsurvivors) and controls (putative survivors) from the most affected population, we carried out a genome-wide association study to identify candidate genes for resistance and susceptibility to CeMV. The genomic data set consisted of 278,147,988 sequence reads and 35,493 high-quality SNPs genotyped across 38 individuals. Association analyses found highly significant differences in allele and genotype frequencies among cases and controls at 65 SNPs, and Random Forests conservatively identified eight as candidates. Annotation of these SNPs identified five candidate genes (MAPK8, FBXW11, INADL, ANK3 and ACOX3) with functions associated with stress, pain and immune responses. Our findings provide the first insights into the genetic basis of host defence to this highly contagious disease, enabling the development of an applied evolutionary framework to monitor CeMV resistance across cetacean species. Biomarkers could now be established to assess potential risk factors associated with these genes in other CeMV-affected cetacean populations and species. These results could also possibly aid in the advancement of vaccines against morbilliviruses.

Phylogenomic history of enigmatic pygmy perches: implications for biogeography, taxonomy and conservation.

Pygmy perches (Percichthyidae) are a group of poorly dispersing freshwater fishes that have a puzzling biogeographic disjunction across southern Australia. Current understanding of pygmy perch phylogenetic relationships suggests past east-west migrations across a vast expanse of now arid habitat in central southern Australia, a region lacking contemporary rivers. Pygmy perches also represent a threatened group with confusing taxonomy and potentially cryptic species diversity. Here, we present the first study of the evolutionary history of pygmy perches based on genome-wide information. Data from 13 991 ddRAD loci and a concatenated sequence of 1 075 734 bp were generated for all currently described and potentially cryptic species. Phylogenetic relationships, biogeographic history and cryptic diversification were inferred using a framework that combines phylogenomics, species delimitation and estimation of divergence times. The genome-wide phylogeny clarified the biogeographic history of pygmy perches, demonstrating multiple east-west events of divergence within the group across the Australian continent. These results also resolved discordance between nuclear and mitochondrial data from a previous study. In addition, we propose three cryptic species within a southwestern species complex. The finding of potentially new species demonstrates that pygmy perches may be even more susceptible to ecological and demographic threats than previously thought. Our results have substantial implications for improving conservation legislation of pygmy perch lineages, especially in southwestern Western Australia.

From conservation genetics to conservation genomics: a genome-wide assessment of blue whales (Balaenoptera musculus) in Australian feeding aggregations.

Genetic datasets of tens of markers have been superseded through next-generation sequencing technology with genome-wide datasets of thousands of markers. Genomic datasets improve our power to detect low population structure and identify adaptive divergence. The increased population-level knowledge can inform the conservation management of endangered species, such as the blue whale (Balaenoptera musculus). In Australia, there are two known feeding aggregations of the pygmy blue whale (B. m. brevicauda) which have shown no evidence of genetic structure based on a small dataset of 10 microsatellites and mtDNA. Here, we develop and implement a high-resolution dataset of 8294 genome-wide filtered single nucleotide polymorphisms, the first of its kind for blue whales. We use these data to assess whether the Australian feeding aggregations constitute one population and to test for the first time whether there is adaptive divergence between the feeding aggregations. We found no evidence of neutral population structure and negligible evidence of adaptive divergence. We propose that individuals likely travel widely between feeding areas and to breeding areas, which would require them to be adapted to a wide range of environmental conditions. This has important implications for their conservation as this blue whale population is likely vulnerable to a range of anthropogenic threats both off Australia and elsewhere.

Genotyping-by-sequencing for estimating relatedness in nonmodel organisms: Avoiding the trap of precise bias.

There has been remarkably little attention to using the high resolution provided by genotyping-by-sequencing (i.e., RADseq and similar methods) for assessing relatedness in wildlife populations. A major hurdle is the genotyping error, especially allelic dropout, often found in this type of data that could lead to downward-biased, yet precise, estimates of relatedness. Here, we assess the applicability of genotyping-by-sequencing for relatedness inferences given its relatively high genotyping error rate. Individuals of known relatedness were simulated under genotyping error, allelic dropout and missing data scenarios based on an empirical ddRAD data set, and their true relatedness was compared to that estimated by seven relatedness estimators. We found that an estimator chosen through such analyses can circumvent the influence of genotyping error, with the estimator of Ritland (Genetics Research, 67, 175) shown to be unaffected by allelic dropout and to be the most accurate when there is genotyping error. We also found that the choice of estimator should not rely solely on the strength of correlation between estimated and true relatedness as a strong correlation does not necessarily mean estimates are close to true relatedness. We also demonstrated how even a large SNP data set with genotyping error (allelic dropout or otherwise) or missing data still performs better than a perfectly genotyped microsatellite data set of tens of markers. The simulation-based approach used here can be easily implemented by others on their own genotyping-by-sequencing data sets to confirm the most appropriate and powerful estimator for their data.

Ecological disturbance influences adaptive divergence despite high gene flow in golden perch (Macquaria ambigua): Implications for management and resilience to climate change.

Populations that are adaptively divergent but maintain high gene flow may have greater resilience to environmental change as gene flow allows the spread of alleles that have already been tested elsewhere. In addition, populations naturally subjected to ecological disturbance may already hold resilience to future environmental change. Confirming this necessitates ecological genomic studies of high dispersal, generalist species. Here we perform one such study on golden perch (Macquaria ambigua) in the Murray-Darling Basin (MDB), Australia, using a genome-wide SNP data set. The MDB spans across arid to wet and temperate to subtropical environments, with low to high ecological disturbance in the form of low to high hydrological variability. We found high gene flow across the basin and three populations with low neutral differentiation. Genotype-environment association analyses detected adaptive divergence predominantly linked to an arid region with highly variable riverine flow, and candidate loci included functions related to fat storage, stress and molecular or tissue repair. The high connectivity of golden perch in the MDB will likely allow locally adaptive traits in its most arid and hydrologically variable environment to spread and be selected in localities that are predicted to become arid and hydrologically variable in future climates. High connectivity in golden perch is likely due to their generalist life history and efforts of fisheries management. Our study adds to growing evidence of adaptation in the face of gene flow and highlights the importance of considering ecological disturbance and adaptive divergence in biodiversity management.

Genome-wide data delimits multiple climate-determined species ranges in a widespread Australian fish, the golden perch (Macquaria ambigua).

Species range limits often fluctuate in space and time in response to variation in environmental factors and to gradual niche evolution due to changes in adaptive traits. We used genome-wide data to investigate evolutionary divergence and species range limits in a generalist and highly dispersive fish species that shows an unusually wide distribution across arid and semi-arid regions of Australia. We generated ddRAD data (18,979 filtered SNPs and 1.725million bp of sequences) for samples from 27 localities spanning the native range of golden perch, Macquaria ambigua (Teleostei; Percichthyidae). Our analytical framework uses population genomics to assess connectivity and population structure using model-based and model-free approaches, phylogenetics to clarify evolutionary relationships, and a coalescent-based Bayesian species delimitation method to assess statistical support of inferred species boundaries. Addressing uncertainties regarding range limits and taxonomy is particularly relevant for this iconic Australian species because of the intensive stocking activities undertaken to support its recreational fishery and its predicted range shifts associated with ongoing climate change. Strong population genomic, phylogenetic, and coalescent species delimitation support was obtained for three separately evolving metapopulation lineages, each lineage should be considered a distinct cryptic species of golden perch. Their range limits match the climate-determined boundaries of main river basins, despite the ability of golden perch to cross drainage divides. We also identified cases suggestive of anthropogenic hybridization between lineages due to stocking of this recreationally important fish, as well as a potential hybrid zone with a temporally stable pattern of admixture. Our work informs on the consequences of aridification in the evolution of aquatic organisms, a topic poorly represented in the literature. It also shows that genome-scale data can substantially improve and rectify inferences about taxonomy, hybridization and conservation management previously proposed by detailed genetic studies.

swinger: a user-friendly computer program to establish captive breeding groups that minimize relatedness without pedigree information.

Captive breeding programmes are often a necessity for the continued persistence of a population or species. They typically have the goal of maintaining genetic diversity and minimizing inbreeding. However, most captive breeding programmes have been based on the assumption that the founding breeders are unrelated and outbred, even though in situ anthropogenic impacts often mean these founders may have high relatedness and substantial inbreeding. In addition, polygamous group-breeding species in captivity often have uncertain pedigrees, making it difficult to select the group composition for subsequent breeding. Molecular-based estimates of relatedness and inbreeding may instead be used to select breeding groups (≥two individuals) that minimize relatedness and filter out inbred individuals. swinger constructs breeding groups based on molecular estimates of relatedness and inbreeding. The number of possible combinations of breeding groups quickly becomes intractable by hand. swinger was designed to overcome this major issue in ex situ conservation biology. The user can specify parameters within swinger to reach breeding solutions that suit the mating system of the target species and available resources. We provide evidence of the efficiency of the software with an empirical example and using simulations. The only data required are a typical molecular marker data set, such as a microsatellite or SNP data set, from which estimates of inbreeding and pairwise relatedness may be obtained. Such molecular data sets are becoming easier to gather from non-model organisms with next-generation sequencing technology. swinger is an open-source software with a user-friendly interface and is available at http://www.molecularecology.flinders.edu.au/molecular-ecology-lab/software/swinger/swinger/ and https://github.com/Yuma248/Swinger.

Towards population-level conservation in the critically endangered Antarctic blue whale: the number and distribution of their populations.

Population-level conservation is required to prevent biodiversity loss within a species, but it first necessitates determining the number and distribution of populations. Many whale populations are still depleted due to 20th century whaling. Whales are one of the most logistically difficult and expensive animals to study because of their mobility, pelagic lifestyle and often remote habitat. We tackle the question of population structure in the Antarctic blue whale (Balaenoptera musculus intermedia) - a critically endangered subspecies and the largest extant animal - by capitalizing on the largest genetic dataset to date for Antarctic blue whales. We found evidence of three populations that are sympatric in the Antarctic feeding grounds and likely occupy separate breeding grounds. Our study adds to knowledge of population structure in the Antarctic blue whale. Future research should invest in locating the breeding grounds and migratory routes of Antarctic blue whales through satellite telemetry to confirm their population structure and allow population-level conservation.

Low genetic diversity in pygmy blue whales is due to climate-induced diversification rather than anthropogenic impacts.

Unusually low genetic diversity can be a warning of an urgent need to mitigate causative anthropogenic activities. However, current low levels of genetic diversity in a population could also be due to natural historical events, including recent evolutionary divergence, or long-term persistence at a small population size. Here, we determine whether the relatively low genetic diversity of pygmy blue whales (Balaenoptera musculus brevicauda) in Australia is due to natural causes or overexploitation. We apply recently developed analytical approaches in the largest genetic dataset ever compiled to study blue whales (297 samples collected after whaling and representing lineages from Australia, Antarctica and Chile). We find that low levels of genetic diversity in Australia are due to a natural founder event from Antarctic blue whales (Balaenoptera musculus intermedia) that occurred around the Last Glacial Maximum, followed by evolutionary divergence. Historical climate change has therefore driven the evolution of blue whales into genetically, phenotypically and behaviourally distinct lineages that will likely be influenced by future climate change.

Hybridization of Southern Hemisphere blue whale subspecies and a sympatric area off Antarctica: impacts of whaling or climate change?

Understanding the degree of genetic exchange between subspecies and populations is vital for the appropriate management of endangered species. Blue whales (Balaenoptera musculus) have two recognized Southern Hemisphere subspecies that show differences in geographic distribution, morphology, vocalizations and genetics. During the austral summer feeding season, the Antarctic blue whale (B. m. intermedia) is found in polar waters and the pygmy blue whale (B. m. brevicauda) in temperate waters. Here, we genetically analyzed samples collected during the feeding season to report on several cases of hybridization between the two recognized blue whale Southern Hemisphere subspecies in a previously unconfirmed sympatric area off Antarctica. This means the pygmy blue whales using waters off Antarctica may migrate and then breed during the austral winter with the Antarctic subspecies. Alternatively, the subspecies may interbreed off Antarctica outside the expected austral winter breeding season. The genetically estimated recent migration rates from the pygmy to Antarctic subspecies were greater than estimates of evolutionary migration rates and previous estimates based on morphology of whaling catches. This discrepancy may be due to differences in the methods or an increase in the proportion of pygmy blue whales off Antarctica within the last four decades. Potential causes for the latter are whaling, anthropogenic climate change or a combination of these and may have led to hybridization between the subspecies. Our findings challenge the current knowledge about the breeding behaviour of the world's largest animal and provide key information that can be incorporated into management and conservation practices for this endangered species.

Microsatellite markers for the praying mantid Ciulfina rentzi (Liturgusidae).

Nine polymorphic microsatellite loci were characterized from an enrichment library of the Australian praying mantid Ciulfina rentzi, a group with a unique reproductive morphology and behaviour. The number of alleles per locus ranged from three to 16 and heterozygosity from 0.24 to 0.94. These markers are the first microsatellites developed for any praying mantid. They will be useful for paternity analysis and for population genetic studies in the Wet Tropics World Heritage Region of Australia.